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Abstract 


This report describes the formulation of a model of the dynamic 
behavior of the Benchmark Active Controls Technology (BACT) wind tunnel 
model for active control design and analysis applications . The model is formed 
by combining the equations of motion for the BACT wind tunnel model with 
actuator models and a model of wind tunnel turbulence. The primary focus of 
this report is the development of the equations of motion from first principles by 
using Lagrange’ s equations and the principle of virtual work. A numerical form 
of the model is generated by making use of parameters obtained from both 
experiment and analysis. Comparisons between experimental and analytical data 
obtained from the numerical model show excellent agreement and suggest that 
simple coefficient-based aerodynamics are sufficient to accurately characterize the 
aeroelastic response of the BACT wind tunnel model. The equations of motion 
developed herein have been used to aid in the design and analysis of a number of 
flutter suppression controllers that have been successfully implemented. 


Introduction 

Active control of aeroelastic phenomena, especially in the transonic speed regime, is a key technology 
for future aircraft design .W The Benchmark Active Controls Technology (BACT) project is part of NASA 
Langley Research Center’s Benchmark Models Program^ 1 * 2 ! for studying transonic aeroelastic phenomena. 
The BACT wind tunnel model was developed to collect high quality unsteady aerodynamic data (pressures 
and loads) at transonic flutter conditions and demonstrate flutter suppression by using spoilers (alone and in 
concert with a traditional trailing edge control surface). The availability of truly multivariable control laws 
also provided an opportunity to evaluate the effectiveness of a controller performance evaluation (CPE) 
tool^ to assess open- and closed-loop stability and controller robustness when applied to multivariable 
systems. An underlying requirement of these objectives was the availability of a mathematical model of the 
BACT dynamics. 

A mathematical model is the basis for nearly all control design methods, therefore an appropriate model 
is essential. The importance of having a good model of the dynamic behavior which satisfies the needs of 
control law design and to accurately assess system performance and robustness cannot be overstated. In 
addition, a model is required to perform the extensive analysis and simulation that are usually required before 
controller implementation to assure that safety is not compromised. This is especially true in the area of 
aeroservoelastic testing in which failure can result in destruction of the wind tunnel model and damage to 
the wind tunnel. 

One of the most important properties to accurately model is the frequency response in the vicinity of 
the key dynamics over the anticipated range of operating conditions. In the case of flutter suppression, the 
key dynamics occur near the flutter frequency and the operating conditions correspond to a wide range of 
dynamic pressures and Mach numbers representing both stable and unstable conditions. The key parametric 
variations associated with the uncertainties due to the assumptions and limitations of the analysis tools and 
other data used to build the model are also important. The development of the mathematical model of the 
dynamic behavior of the BACT presented herein was motivated by these factors. 

Starting from first principles and by using appropriate idealizations of the structural and aerodynamic 
characteristics, a model was developed that includes spoiler controls and explicitly contains the key physical 
parameters. This report emphasizes the development of the equations of motion of the BACT wind tunnel 
model. A more comprehensive mathematical model that includes actuator and turbulence models is presented 
as well. The resulting numerical model is also validated by using several types of static and dynamic 
analysis. 



Nomenclature 


Symbols 


Sq , Z?j , #2 control effectiveness matrices associated with control deflection, rate, and 
acceleration, respectively 

b wing semispan (distance from root to tip) 

C L , C M aerodynamic lift and pitching moment coefficients, respectively 

C L( .), derivative of lift and moment coefficients with respect to. (•) , respectively 


c 

c 

A, 


d 


(•) 


E 




e 2 


e h 


e 6 


g 

H S 

h 

h 

h 

K e , K 6 


wing chord 

mean aerodynamic chord 

structural and aerodynamic damping matrices, respectively 
chordwise distance from the shear center to the (•) sensor location 
wind tunnel turbulence disturbance matrix 

distance from origin of body fixed reference to center of gravity of B ACT/PAPA 
body (excluding control surfaces) 

distance from origin of body fixed reference to center of gravity of trailing edge 
control surface 

distance from origin of body fixed reference to hinge line of trailing edge control 
surface 

distance from hinge line of trailing edge control surface to center of gravity of 
trailing edge control surface 

gravitational acceleration 

aerodynamic control surface hinge moment 

plunge or vertical deflection of wing (positive down) 

angular moment of interia of BACT/PAPA body (excluding control surfaces) 
about the origin of the body fixed reference 

angular moment of inertia of trailing edge control surface about its hinge line 
plunge, pitch, and control surface support stiffness, respectively 


K s , K a structural and aerodynamic stiffness matrices, respectively 

L aerodynamic lift force 

Lj turbulence model parameter 

i{x) chordwise distance form the origin of the body fixed coordinate system to the 

point at which angle of attack is referenced 

M aerodynamic pitching moment, Mach number 
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M s , M a , M g structural and aerodynamic mass matrices, respectively 


W2, 

m 2 

p(x,y,t ) 

Q() e 

Q() nc 

Q o 
Qt 
q 
q 

if 

s 

s ij 

T 


u e 

U 0 

sw 


W 

JC 

x 

X 


p 

f 


y 

z 

a 

h 

y P 


mass of B ACT/PAPA body (excluding control surfaces) 
mass of trailing edge control surface 

differential pressure (upper surface - lower surface) at the point (x ,y) on the wing 
surface at time t 

generalized externally applied force associated with generalized coordinate (•) 

generalized nonconservative (i.e M structural damping) forces 

generalized aerodynamic force matrix at zero angle of attack 

generalized aerodynamic force matrix associated with turntable angle 

pitch rate 

dynamic pressure 

flutter dynamic pressure 

planform area of wing 

inertial coupling between generalized coordinates i and j 

kinetic energy 
time 

gravitational potential energy 
elastic strain energy 
freestream airspeed 

work done by virtual displacements of generalized coordinates 

normal perturbation velocity of local flow field disturbances (positive down 
relative to the freestream flow) 

vector of disturbances 

chordwise distance (positive aft from root mid-chord) 
chordwise position of origin of body fixed coordinate system 

chordwise position of hinge of control surface 

spanwise distance (positive outboard from wing root) 
displacement normal to wing chord (positive up) 
angle of attack, turbulence model parameter 
turbulence model parameter 

turbulence model parameter 
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6 control surface deflection, vector of control surface deflections 


6 C 

commanded control surface deflection 

&TE’ &US 

trailing edge and upper spoiler control surface deflections, respectively 

Vg 

random input to turbulence model 

S*. Cg 

modal structural damping ratio, plunge and pitch, respectively 

e 

pitch angle 

O'f 

turntable angle 

«<■>• § 

generalized coordinate, vector of generalized coordinates 

(O h , co 0 

in vacuo modal frequencies, plunge and pitch, respectively 

COf 

flutter frequency 

Subscripts / Superscripts 

Oo 

equilibrium, trim, or bias value of (*) 

( : ) 

perturbation of (*) 

(•)/ 

flutter 

()l£/ 

leading edge inboard 

(')leo 

leading edge outboard 

(')re 

trailing edge control 

(')tei 

trailing edge inboard 

(')teo 

trailing edge outboard 

(■)t/S 

upper spoiler control 

(Y 

external 

(■r 

nonconservative 

Operators 


(0 

derivative with respect to time 

(•) 

second derivative with respect to time 

d 

partial derivative with respect to. (•) 

Abbreviations 


BACT 

Benchmark Active Control Technology 
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CFD 


Computational Fluid Dynamics 
CPE Controller Performance Evaluation 

FFT Fast Fourier Transform 

NACA National Advisory Council for Aeronautics 

PAPA Pitch and Plunge Apparatus 

psf pounds per square foot 

rms root mean square 

TDT Transonic Dynamics Tunnel 

2-DOF two degree of freedom 

The BACT Wind Tunnel Model 

The BACT wind tunnel model is a rigid rectangular wing with an NACA 0012 airfoil section. 14 ^ A 
photograph of the model with physical dimensions is shown in figure 1 . The wing is equipped with a 
trailing edge control surface and upper and lower surface spoilers that can be controlled independently via 
hydraulic actuators. The lower surface spoiler (not shown) is identical in size and relative location to the 
upper surface spoiler. The wind tunnel model is instrumented with pressure transducers, accelerometers, 
control surface position sensors, and hydraulic pressure transducers. The accelerometers are the primary 
sensors for feedback control and are located at each comer of the wing. 



Figure 1 - BACT wing section with dimensions. 


5 



The wing is mounted to a device called the Pitch and Plunge Apparatus^ (PAPA) which is designed to 
permit motion in principally two modes — rotation (or pitch), and vertical translation (or plunge). A 
photograph of the BACT wing mounted to the PAPA is shown in figure 2. The mass, inertia, and center of 
gravity location of the system can be controlled by locating masses at various points along the mounting 
bracket. The stiffness properties can be controlled by changing the properties of the support rods. The 
PAPA is instrumented with strain gauges to measure normal force and pitching moment and is mounted to 
a turntable that can be rotated to control the angle of attack of the wing. 

The combination of the BACT wing section and PAPA mount will be referred to as the BACT system. 
The BACT system was precisely tuned to flutter within the operating range of the NASA Langley 
Transonic Dynamics Tunnel (TDT)^ in which the system was tested. The range of Mach numbers and 
dynamic pressures over which flutter occurs permits the study of transonic aeroelastic phenomena. More 
detailed descriptions of both the BACT wing section and the PAPA mounting system can be found in 
references 4, 6, and 7. 



Figure 2 - Photograph of BACT wing section and PAPA mount. 


Aeroelastic Equations of Motion 

The BACT system has dynamic behavior very similar to the classical two degree of freedom (2-DOF) 
problem in aeroelasticity This similarity was exploited in the development of the aeroelastic equations of 
motion for the BACT system by representing it as a 2-DOF system and by using a strip theory 
aerodynamic approximation J 9,10 ^ The difference between the classical 2-DOF system and the BACT system 
is primarily the complexity of aerodynamic behavior and presence of additional structural modes. The finite 
span and low aspect ratio of the BACT wing introduce significant three dimensional flow effects. The finite 
span of the control surfaces and their close proximity also introduce significant aerodynamic effects. Higher 
frequency structural degrees of freedom are associated with the PAPA mount and the fact that the wing 
section is not truly rigid. These effects are accounted for by various means and will be discussed 
subsequently . 

Lagrange’s equations were used to derive the equations of motion for the BACT system and the 
principle of virtual work was used to obtain expressions for viscous damping and the generalized 
aerodynamic forces. Lagrange’s equations and the principle of virtual work provide a simple and straight 
forward method for deriving the equations of motion for aeroelastic systems Lagrange’s equations readily 
allow one to represent motions relative to a moving frame. The principle of virtual work has the advantage 
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of automatically accounting for the forces of constraint and thereby greatly simplifying the determination of 
generalized forces. 

The basic requirement for applying Lagrange’s equations is that the velocity of each point in the body 
be represented in an inertial reference frame. The efficient application of Lagrange’s equations can be 
facilitated by representing the inertial quantities in convenient coordinate systems and selecting an 
appropriate set of generalized coordinates. Once the generalized coordinates are chosen, expressions for the 
kinetic and potential energies are determined and Lagrange's equations are applied. The principal of virtual 
work is used to compute the generalized forces acting on the body. The complete set of equations of motion 
are then assembled. Finally, the equations of motion are used to determine the equilibrium solution and 
perturbation equations for the system and a numerical model appropriate for control system design is 
obtained. 

Coordinate Axes and Generalized Coordinates 

The BACT system can be idealized as a collection of four rigid bodies corresponding to each of the 
three control surfaces and the remaining wing/PAPA element. Figure 3 depicts the relevant quantities for 
the wing and the trailing edge control surface. The spoiler control surfaces were treated in an analogous 
fashion but are omitted here for ease of discussion. 



Figure 3 - Structural representation showing trailing edge control only. 


There were five coordinate systems used. One coordinate system was fixed to the wing and moves with 
it. Another coordinate system is fixed in inertial space and oriented relative to the turntable to which the 
BACT system was mounted. The inertial coordinate system was chosen to coincide with the undeformed 
position of the body-fixed system. The three other coordinate systems are fixed to each of the control 
surfaces and rotate relative to the body-fixed system about the hinge of each surface. 

The origin of the inertial coordinate axes is located at the shear center of the undeformed position. The 
origin of the body fixed coordinate axes coincides with the instantaneous shear center of the system. The 
origin of the control surface-fixed coordinate axes coincides with the hinge lines. 

The generalized coordinates were selected to simplify the derivation of the equations of motion and were 
based on some key assumptions about the nature of the motion of the BACT system. The wing section and 
each control surface are assumed to be rigid in both the spanwise and chordwise directions. It is also 
assumed that the wing motion is limited to two degrees of freedom — pitch and plunge. This assumption 
implies that the other structural modes of the BACT system are insignificant. Investigation of the structural 
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vibration characteristics of the PAPA mount with a very similar wing model was shown to support this 
assumption. 11 ^ The next lowest frequency for any transverse mode was more than six times the frequency 
of the pitch and plunge modes and outside the frequency range of interest. 

Based on these assumptions the BACT requires five generalized coordinates, i=I,2,...,5. The two 
associated with the pitch and plunge degrees of freedom were chosen to be pitch angle 6 and plunge 
displacement h of the body fixed coordinate axes relative to the inertial coordinate axes. The three associated 
with the angular rotation of the control surfaces were chosen to be the trailing edge, upper spoiler, and 
lower spoiler control surface angles (5 -ff , ^US » ^ndd^, respectively. The three coordinates h, 0, and 
<5 * 6 te are depicted in figure 3 for the system in which only the trailing edge control surface is included. 

Kinetic and Potential Energy 

The selection of the generalized coordinates allows expressions for the kinetic and potential energies to 
be formulated. The kinetic energy of a body is the work required to increase its velocity from rest to some 
value relative to an inertial frame. By using the quantities defined in figure 3 the kinetic energy expression 
for the BACT system neglecting the spoiler controls can be written 

T = — m^ih + e^d) 2 + — l\6 2 + — m 2 {h + e 2 d + + T ^2^ + • 0) 

2 2 2 2 

The potential energy of a rigid body supported by pitch and plunge springs consists of two terms, the 
strain energy in the springs and the gravitational potential energy. The gravitational potential is defined 
relative to a arbitrary datum. By using the quantities defined in figure 3 and assuming the datum to be the 
origin of the inertial frame the gravitational potential for the BACT system is 

U g = + e\ sin 0 )cos 07 - - m 2 g(h + e 2 sind + e 6 sind)cosd T , (2) 

where g is the gravitational acceleration and d T is the turntable angle. 

The strain energy is the work done in going from the undeformed reference position to the deformed 
position. By using the quantities defined in figure 3 the strain energy for the BACT system can be written 

U e = X - K h h 2 + ~ K e d 2 + i K S {6 - 6 C ) 2 , (3) 

where K h and Kg are the spring constants associated with the stiffness properties of the PAPA mount and 
K s is the spring constant representing the flexibility in the structure supporting the actuator. Note that the 
strain energy associated with the control surface is chosen to be based on the difference between the actual 
and the commanded control surface rotations <5 and 6 C , respectively. 

The total potential energy for the BACT system is simply the sum of the gravitational potential and 
the strain energies in Eqns (2) and (3), respectively. 

Applying Lagrange’s Equations 

Lagrange’s equations can be expressed as 

d ( jt\ sr su _ ,,, 
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where T and U are the total kinetic and potential energies of the system, is the ith generalized 
coordinate, and Q ^ is the generalized force associated with ^ and includes externally applied forces, 

nonconservative forces, and forces of constraint. 

Applying the expressions for the kinetic and potential energies to Eqn (4) results in the system of 
equations 
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s h6 COS <5 


Qs 


(5) 


The terms m , Iq , and are the generalized masses of the pitch, plunge, and control surface modes, 
respectively. The terms and are the inertial coupling between the generalized coordinates. 

These terms are related to the quantities defined in figure 3 by 

m s wi] + m2 , 

Iq = 7 | + I2 + ^ + ftl'^6'2^ , 

I 6 =l 2 +m 2 e2 2 , (6) 

She sm l e l + " I 2 e 2> 

s h6= m 2 e 6> 
s 9d = 1 2 + m 2 e 2 e 6 . 


If one assumes that the control surface stiffness is very large (i.e., the deformation due to hinge load is 
insignificant) then Eqn (5) can be simplified by eliminating the generalized force associated with the control 
surface Q$ . Assume that the control surface stiffness is very large so that 


* -6 


\_ 
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e « 1 . 


After eliminating terms of order e, Eqn (5) can be approximated by 
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( 8 ) 


Note that the value of 6 has been assumed to be small so that cos0 is approximately unity. Also note that 
the inertial coupling between the wing structure and the control surface is retained in the equations. 

All that remains to complete the equations of motion is to determine expressions for the generalized 
forces Qh and Qq . 

Applying the Principle of Virtual Work 

The principle of virtual work can be applied to obtain expressions for the generalized forces. The basic 
advantage of using this method is that the forces of constraint are eliminated automatically. In addition, the 
principle of virtual work can be used to determine expressions for dissipative forces such as damping. 
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The generalized force Q ^ can be determined from 


<k- 


dSW 


( 9 ) 


where 6W is the work done on the system by arbitrary infinitesimal (or virtual) displacements of the 
generalized coordinates <5£. This virtual work includes the work done by nonconservative forces (e.g., 
damping) and external forces. The work done by forces of constraint are zero under virtual displacements. 

Nonconservative ( Damping ) Forces 

Structural damping is often characterized as a viscous force. Experimental data suggests that this is a 
reasonable assumption for the BACT system undergoing small motions. Viscous forces are those where 
the force varies in proportion to the velocity at the point the force is applied but in the opposite direction. 
The generalized damping forces for the BACT can be represented as 


Qh nc ] 
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The terms in Eqn (10) were chosen so that the damping coefficients associated with the plunge and pitch 
modes, and respectively, correspond to those obtained from experiment. The matrix premultiplying 
the diagonal damping matrix is the mass matrix from Eqn (8). The other constants a) h and 0)q correspond 
to the in vacuo vibration frequencies for the pitch and plunge modes. This so-called equivalent modal 
damping formulation results in a diagonal damping matrix when the system is written in modal form. 

External ( Aerodynamic ) Forces 

The externally applied forces are due to aerodynamics and result from the distributed pressures applied to 
the surface of the BACT wing. The virtual work for the wing can be written 

<5VV = J^J*^p(x,y,t)dz(x,y,t) dxdy , (11) 

where b is the wing semi-span, c is the section chord, p(x,y,/) is the differential pressure distribution over 
the surface of the wing, and dz(x,yj) is a virtual displacement normal to the wing surface. Figure 4 
depicts the pressure distribution for a chordwise section. 



Figure 4 - Pressure distribution over wing section. 
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The virtual displacement 6z can be written in terms of virtual displacements of the generalized 
coordinates h, 6 , and 6 as 


«>+<*-V*fl(t)+j (j[ _^ )&5 , Wx,- <12) 

Substituting this expression in Eqn (11) and performing the differentiation described in Eqn (9) results in 
expressions for the generalized aerodynamic forces 

Qh = ^ = ~sif. c ^p(x>y,t)dxdy s -L, 

Qe e *^T = _ Ja £^P(x,y,t)(x-x p )dxdy = M p , (13) 

Qb = ^ = -fifxpPit’y^ ~x f )dxdy = H 6 , 


where L, M p , and are the lift, pitching moment about the reference point x p and the control surface 
moment about its hinge line x j , respectively. The reference point x p was chosen to correspond to the 
shear center (i.e., the origin of the body fixed coordinate system). 

There are several ways in which the integrated surface pressures can be approximated in practice. If 
computational aerodynamic analysis results are available the integration can be approximated by using the 
pressures on the computational grid. If experimental force and moment data are available it can be used 
directly in Eqn (13) after accounting for differences in the moment reference point. 

A common method of approximating the aerodynamic forces is to use stability and control derivatives. 
The aerodynamic forces are represented as a linear function of angle of attack, control surface deflection, and 
their rates by 


L = qSC L 


= qS 


C Lq + C L a a + C L r & + 


2Ur 


-(c La a + c Lq 6 + c L& 6) 


(14) 


M = qScC M 


= qSc 


C M 0 + C M a a + C M 6 & + 


2Ui 


■ ( ® + C M q 0 + 


(15) 


A description of each coefficient is presented in table 2 in the Numerical Model Parameters section. A 
comparable expression for the hinge moment has been omitted since the need for the generalized force 
associated with hinge moment in the equations of motion was eliminated in Eqn (8). 

The coefficient based approach is a rather simplistic way to represent the aerodynamic forces and 
moments, but it is quite acceptable for the BACT system as will be seen. The need to include more 
sophisticated aerodynamic modeling approaches such as rational function approximations or time accurate 
CFD is mitigated by the fact that the reduced frequency for the BACT system is relatively low (e.g., 
approximately 0.04 at a flutter dynamic pressure of 150 psf using semichord as the reference length). 
Therefore, the lag between changes in angle of attack and the resulting lift force is not significant. 
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A relation between the generalized coordinates and the angle of attack is required to represent the 
generalized aerodynamic forces in a form consistent with the equations of motion. Here it is assumed that 
the wind tunnel model is rigid and that the angle of attack is the same for each spanwise position. Based on 
the choice of the generalized coordinates, the angle-of-attack at any chordwise point x at any time t can be 
represented by 


, ,, „ a . , h(t ) ^(Jt)0(O 

a(x,t) = 9 t + d(t) + + — - - ^ , 

u 0 u 0 u 0 


(16) 


where #7 is the turntable angle, l(x) is the chordwise distance from the origin of body fixed coordinate 
system to the point at which the angle of attack is referenced (positive aft), w g is the normal perturbation 
velocity of disturbances in the local flow field (positive down relative to the free stream flow), and ( 7 0 is 
the freestream velocity. 

By using the expression for angle of attack from Eqn (16) in Eqns (14) and (15) the generalized applied 
external forces can be written 
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(17) 


Complete BACT Equations of Motion 

Combining the generalized forces, Eqns (10) and (17), with Eqn (8) results in a complete set of 
equations of motion for the BACT system. The general form of the equations can be expressed as 

( Ms ~2U^ M ^f + [ Ds ~^ Da )% + ( K s-4 SK °^ 

_ _ ( 18 ) 

~ qSQ§ e + qSQjQj + M g g cosO^ + f?2<5 + — — B\8 + qSBQ& + - — Ew , 

2U 0 U 0 

where § is the vector of generalized coordinates (see Appendix), 0 T is the turntable angle, <5 is the vector 
of control surface deflections, w is the vector of disturbance inputs (downwash velocity and acceleration - 
see Appendix), and q is the dynamic pressure. The definitions of the matrices M s , M a , M gf D s , D a , 

K a , K a> Q 0j Q t , B 0 , Z?j, B 2 ,and E are presented in the Appendix. 

Notice that the effect of the aerodynamic forces is to modify the mass, damping, and stiffness properties 
of the system. It is this aerodynamic coupling that is the essential feature of aeroelastic systems and leads to 
the flutter instability. Also note that the aerodynamic matrices are explicit functions of the dynamic 
pressure and freestream velocity. So, for a given Mach number, the equations of motion represent a 
continuum of system dynamics associated with a range of dynamic pressures. This parameterization is 
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particularly useful for studying the stability characteristics of the system. Finally, note that there are three 
terms in the equations that are constant, assuming the turntable angle is fixed. These terms determine the 
static equilibrium of the system. 

Equilibrium Solution and Perturbation Equations 

The static equilibrium of the BACT system is obtained by setting all time derivatives in Eqn (18) to 
zero and by solving for the generalized coordinates. Doing so results in 

%0={ K s- K a) \Qo € +Qt6t + M g g COS #7* + V»o]- (I 9 ) 


The subscript (*)o on the generalized coordinate vector is used to indicate static equilibrium. The subscript 
(*)o ° n the control input is used to denote a fixed reference value (e.g., bias, trim value). 

The generalized coordinates can be expressed as the sum of the static (or equilibrium) part a 

perturbation part £ by 


§-§0+t- ( 2 °) 

The control input can be expressed as the sum of the bias or static part and the time varying or dynamic 
part 6 by 


<5 = <5q + 3 


(21) 


Substituting Eqns (20) and (21) into Eqn (18), by using the fact that £q, 6 ^, and 60 are constant, and 
by eliminating the constant terms by using Eqn (19) results in the perturbation state space equations of 
motion for the BACT system 
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2Ui, 


M a r‘ 1—E 

Un 


( 22 ) 


While the form of the equations that appear in Eqn (18) describe the complete motion of the system, it 
is the form of the equations of motion presented in Eqn (22) that are most readily applicable to typical 
control system design methods. Note that even though Eqn (22) was derived for a single control surface, 
extension to multiple control surfaces is straight forward since there is no inertial coupling between the 
various control surfaces of the BACT. There is, however, aerodynamic coupling between the control 
surfaces due to their close proximity. The form of the aerodynamic force expressions allows some effects of 
coupling to be approximated within the stability and control derivative terms by altering the derivative 
values to account for control surface biases. For example, the effectiveness of the trailing edge control 
surface ( C Lte , C ) can be adjusted to account for the influence of biases of the other surfaces. However, 

perturbation effects (i.e., changes in the control effectiveness due to motion of the other surfaces about their 
bias values) are ignored in this approach. 
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Numerical Model Parameters 


A numerical form of the equations of motion of the BACT system was obtained by substituting 
numerical values for each parameter in the equations of motion developed above. Most of the parameter 
values were obtained from experimental data but some of the aerodynamic data were obtained from 
numerical analysis using computational aerodynamics. 

The mass and inertia parameters were obtained by measuring the mass, stiffness, and damping 
properties of the various components of the BACT system. The geometric parameters (e.g., centers of 
gravity, shear center, sensor locations, and aerodynamic reference quantities) were also obtained directly from 
measurement of the BACT and PAPA components. The mass, stiffness, and damping parameter values that 
were used in the numerical form of the BACT equations of motion are presented in table 1 . The center of 
gravity and shear center were essentially coincident and located at the mid-chord point of the wing. 

The static aerodynamic parameters were determined from experimental data that were obtained from a 
previous wind tunnel test in which the BACT wing was mounted on a force and moment balance J 7 ^ Force 
and moment data for various angles of attack and control surface positions were used to compute most of 
the stability and control derivatives using finite differences. The upper spoiler control derivatives were 
obtained relative to a bias deflection (5 to 10 degrees) about which the spoiler was deflected to control the 
system. The spoiler bias allows both positive and negative control inputs to be produced. 


Table 1 - Mass, Stiffness, and Damping Parameters 


Symbol 

Description 

Value 

Units 

m 

mass 

6.08 

slug 

h 

pitch inertia 

2.80 

slug-ft 2 

Kh 

plunge stiffness 

2686 

lb/ft 

K e 

pitch stiffness 

3000 

lb-ft 

0) h 

in vacuo plunge frequency 

21.01 

rad/sec 

«>e 

in vacuo pitch frequency 

32.72 

rad/sec 


plunge damping ratio 

0.0014 

- 

Ze 

pitch damping ratio 

0.0010 

- 

s hd 

pitch-plunge inertial coupling 

0.0142 

slug-ft 

s h&TE 

plunge-TE inertial coupling 

0.00288 

slug-ft 

s B5te 

pitch-TE inertial coupling 

0.00157 

slug-ft 2 

Sh &us 

plunge-US inertial coupling 

0.00039 

slug-ft 

S 86us 

pitch-US inertial coupling 

9.8e-05 

slug-ft 2 
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Sufficient experimental data were only available to quantify the trailing edge and upper spoiler control 
surface aerodynamic characteristics. Available data for the lower spoiler were not complete enough to 
characterize its aerodynamics. In addition, there were insufficient data to account for aerodynamic coupling 
between the spoiler and trailing edge control. Therefore, the numerical model is limited to the trailing edge 
and upper spoiler surfaces with no aerodynamic coupling between controls (i.e., the influence of spoiler bias 
on trailing edge control effectiveness is ignored). 

The dynamic derivatives (e.g., C M < and C L ^ ) were obtained from computational aerodynamic 

analysis. Analytical values were used because they were available from models previously generated using 
ISAC.* However, the dynamic derivatives associated with spoiler-surface deflection rate and acceleration 
were not available from these models and so all the control surface rate and acceleration derivatives were 
assumed to be zero. 

The numerical values for the static and dynamic stability and control derivatives are presented in 
table 2. These values were determined at a single Mach number of 0.77 and a single dynamic pressure of 
143 psf. The moment coefficients are referenced relative to the shear center that coincides with the mid-chord 


Table 2 - Aerodynamic Parameters 


Symbol 

Description 

Value 

Source 

C *o 

lift at zero angle of attack 

0 

experiment 

C M 0 

pitching moment at zero angle of attack 

0 

experiment 


lift curve slope 

4.584 

experiment 

C M a 

moment curve slope 

1.490 

experiment 


plunge damping due to angle of attack rate 

-3.1054 

analysis 


plunge damping due to pitch rate 

2.5625 

analysis 

C M a 

pitch damping due to angle of attack rate 

-2.6505 

analysis 

C M q 

pitch damping due to pitch rate 

-0.4035 

analysis 

C l& te 

TE lift effectiveness 

0.63 

experiment 

^ M Ste 

TE moment effectiveness 

-0.0246 

experiment 

Cl >us 

US lift effectiveness 

0.22 

experiment 


US moment effectiveness 

0.0573 

experiment 

£ 

distance between shear eerier and aerodynamic center 

-0.175 c 

experiment 

S 

planform area, (ft 2 ) 

3.55 

experiment 

C 

mean aerodynamic chord, (ft) 

1.33 

experiment 


* ISAC (or Integration of Structures, Aerodynamics, and Controls) is a tool for modeling aeroelastic 
aircraft for application to control system analysis and design J 12 ^ 
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point of the wing. Finally, note that the expression for the generalized aerodynamic forces in Eqn (17) 
requires the selection of an angle of attack reference point ^(jt) (i.e., the distance between the shear center 
and the point at which the angle of attack of the wing is measured). The aerodynamic center location based 
on the experimental lift and moment curve slopes was used as the angle of attack reference. This point 
varied with Mach number but was near the 30% chord location at Mach 0.77 and for most other operating 
conditions as well. 

The numerical form of the matrices that appear in Eqn (18) that use the values presented in tables 1 
and 2 are presented in the Appendix. Note that since the model is parameterized by q and Uq , a continuum 
of numerical state space models of the form in Eqn (22) can be generated. However, these models are all 
associated with the single Mach number at which the aerodynamic data were obtained. 


Actuator and Turbulence Models and the Output Equation 

The equations of motion alone are not sufficient to describe the dynamic behavior of the BACT system. 
While the 2-DOF system structure is sufficient to describe the basic aeroelastic properties of the BACT 
system, additional elements are necessary to develop a model suitable for control system design. The 
relative magnitude of the dynamic response is determined by the nature of the disturbance environment and 
influences the control activity required to achieve the desired level of closed-loop performance. Therefore, a 
characterization of the turbulence environment in the TDT is needed, i.e., a turbulence model. The ability of 
the control surfaces to produce the desired activity is dependent on the dynamic response characteristics of 
the actuators including bandwidth, position and rate limits, and other nonlinearities. Therefore, 
characterizations of the actuator dynamics are also needed, i.e., actuator models. Finally, a set of 
measurement signals is required to provide the basis for feedback control. An output equation relating the 
generalized coordinates to the measurement variables is therefore required. It is the combination of the 
actuator models, the turbulence model, output equation, and the aeroelastic equations of motion that 
determine the degree to which a control system will be able to achieve a desired level of performance and 
robustness. 

Actuator Models 

Actuator models of the BACT wind tunnel model were obtained from experimental data by using a 
simple parameter estimation process described in reference 13. This process was used to select the 
parameters of a second order actuator model to minimize the frequency response error over the frequency 
range of interest. The actuator transfer function model is 

<5(s) ka) 2 

TTT = 7 > (23) 

S 2 +2£ a (D a S + Q) a 2 

where k is a gain, and o> a and £ a are frequency and damping, respectively. The parameter values resulting 
from the parameter estimation process are presented in table 3. A state space form of the second order 
actuator model and the numerical matrices for the trailing edge control surface and the upper spoiler 
actuators are presented in the Appendix. 

Other factors that have a significant affect on control system design are the position and rate limits of 
the actuators and nonlinearities such as dead zone and backlash. The position limits for the trailing edge 
control surface and the spoiler-surfaces were approximately 12 degrees and 45 degrees , respectively. The rate 
limits, dead zone, and backlash have not yet been determined. 
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Table 3 - BACT Actuator Model Parameters 


Symbol 

Description 

Value 

Units 

TF. 

TE actuator gain 

1.02 

deg/deg 

t’TE 

TE damping ratio 

0.56 

- 

0J TE 

TE frequency 

165.3 

rad/sec 

k US 

US actuator gain 

1.16 

deg/deg 

£us 

US damping ratio 

0.85 

- 

, £&£ . 

US frequency 

164.0 

rad/sec 


Turbulence Model 


A model of the turbulence environment within the TDT was developed in reference 13 using power 
spectrum data. The model structure is that of a Dryden spectrum with the parameters adjusted to 
approximate the desired power spectral density. The transfer function form of the turbulence model is 


( 2 Jt\ 

Wg(s) _ 2 ) 

Y ( 2 4 jt 4 Jt 2 \ 

I s 

where 



The output w g is the down wash disturbance velocity that appears in Eqns (16) and (17). The input t] g 

is a normally distributed univariate random variable. Table 4 presents a range of values for the turbulence 
model parameters. Note that the parameter values are based on data collected in an air medium, not the 
medium in which the BACT was tested (i.e., R-12) and so the reference speed, Uq, corresponds to a 
different Mach number. No data at the appropriate operating conditions is available. The preferred reference 
speed value is 400 fps because it is close to the test conditions and correlates best with experimental data. 
The preferred value was used to perform the disturbance analysis presented in the next section. A state space 
form of the turbulence model and the corresponding numerical matrices are presented in the Appendix 


Table 4 - TDT Turbulence Model Parameters 



Reference Speed, U 0 (fps) 

Parameter 

100 

200 

300 

400* 

a 

0.01 

0.025 

0.007 

0.082 

h 

0.477 

0.475 

0.521 

0.667 

Yp 

0.546 

0.464 

0.497 

0.533 

—h— 

3.261 

3.71 

3.391 

4.163 


Preferred reference speed value for analysis and simulation. 
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Output Equation 


The BACT system has four accelerometers, one mounted in each comer of the rectangular wing. These 
accelerometers sense vertical acceleration measured in g’s, positive up (opposite to the sign convention for 
plunge, h). The acceleration at any point on the BACT wing, excluding control surfaces, has the form 


z(x,y) = 



(25) 


where h and 6 are the generalized coordinates, d i is the chordwise distance from the shear center to the 
location of the ith sensor, and g is the gravitational acceleration. Table 5 presents the chordwise distance 
(positive aft) between the shear center and the accelerometer for each of the four accelerometers — leading 
edge inboard (LEI), leading edge outboard (LEO), trailing edge inboard (TEI), and trailing edge outboard 
(TEO). 


Table 5 - Accelerometer Locations 


1 

Accelerometer Location i 

.. | 

^LEJ 

d LEO 

d TEI 

d TEO 

Distance (ft) 

i 

-0.599 

-0.599 

0.433 

0.420 


All the components described above were combined to form the complete numerical model of the 
BACT system. The following section addresses a variety of analyses that were performed to assess the 
accuracy and validity of the numerical model . 


Validation of Numerical Model 

There are many ways to assess the validity of the BACT numerical model. In this section a few 
comparisons are made between the properties of the numerical model and the actual BACT wind tunnel data. 
These assessments can be broken down into two categories — static properties and dynamic properties. The 
static properties assess the characteristics of the equilibrium solutions. The dynamic properties assess the 
key response characteristics of the system in the context of flutter behavior and frequency response. 

Static Properties 

The equilibrium position (pitch and plunge) of the BACT system depends on the turntable angle arri 
wind tunnel operating conditions and represents a balance between the elastic and aerodynamic forces acting 
on the wing. Good agreement between the equilibrium position of the wind tunnel model and the 
equilibrium solution of the numerical model would indicate that the stiffness (structural and aerodynamic) 
and control surface effectiveness properties are well modeled. 

Parameters that were recorded during the wind tunnel tests include turntable angle, control surface 
biases, and pitch angle (plunge position was not measured). In addition, the test conditions (i.e., Mach 
number, dynamic pressure, and control surface positions) were recorded. By using these quantities in 
Eqn (19), one can determine the equilibrium pitch and plunge position of the BACT system for comparison 
with experimental data. 

Table 6 presents the computed and measured pitch angle for a small representative set of test 
conditions. The error in equilibrium pitch angle is less than 5 percent for all but one point. In addition, the 
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trends are consistent; increasing the turntable angle increases the pitch angle, deflecting the upper spoiler 
upward decreases the pitch angle, and deflecting the trailing edge control downward decreases the pitch angle. 

Table 7 presents the computed and measured pitch angle increments due to the changes in turntable 
angle and control surface deflections from table 6. The size of the pitch angle increments due to turntable 
angle and spoiler deflection are consistent. The main discrepancy is that the pitch angle increment due to 
trailing edge control deflection from the experimental data is five times higher than from the numerical 
model. The source of the discrepancy could not be identified due to a lack of data. However, the difference 
could be due to experimental error or inaccurate modeling of trailing edge control effectiveness and/or 
stiffness of the PAPA mount (the parameters involved in balancing the aerodynamic and structural pitching 
moments). Experimental error is a more likely cause for the differences because the trailing edge control 
derivatives and the structural stiffness that were used in the numerical model were determined from 
experimental data and should accurately characterize the static properties of the system. 


Table 6 - Static Equilibrium Position Comparison 


Mach 

Number 

Dynamic Pressure 
(psf) 

| e T 
j (deg) 

8 te 

(deg) 

6 US 

(deg) 

B exp 

(deg) 

B model 
(deg) 

Ad 

(%) 

.65 

112 

| 1.6 

0 

0 

2.1 

2.17 

3.3 

.65 

115 

1.6 

0 

-10 

2.0 

2.05 

2.5 

.70 

126 

1 1.6 

\ 

0 

0 

2.4 

2.28 

-5.0 

.70 

126 

1.6 

10 

0 

2.0 

2.20 

10.0 

.77 

120 

1.4 

0 

0 

2.0 

1.95 

-2.5 

.77 

120 

4.5 

0 

0 

6.0 

6.27 

4.5 


Table 7 - Static Equilibrium Sensitivity Comparison 


Mach 

Number 

Dynamic 
Pressure (psf) 

A Q'p 
(deg) 

A8 T e 

(deg) 

A& us 

(deg) 

AB exp 

(deg) 

model 

(deg) 

.65 

112 

0 

0 

10 

0.1 

0.12 

.70 

126 

0 

-10 

0 

0.4 

0.08 

.77 

120 

-3.1 

0 

0 

-4.0 

-4.32 
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Dynamic Properties 


A major concern in using the numerical model for control system design is the ability of the model to 
accurately represent the transition to flutter including the frequency and dynamic pressure of the flutter 
onset. Another concern is the fidelity of the model from a frequency response perspective because of the 
relationship between the frequency response and the design of the control system. Finally, the level of 
response due to turbulence is an important issue because the maximum level of control activity depends 
directly on the level of response of the wind tunnel model to turbulence. Each of these properties were 
reviewed and compared with experimental data (and numerical models generated with IS AC where possible) 
to assess the validity of the numerical model. Note that the IS AC-based models were not corrected using the 
static structural and aerodynamic data that was used to obtain the numerical model discussed herein. 

Flutter Properties 

The BACT wind tunnel model experienced flutter in the TDT at a dynamic pressure of approximately 
148 psf at a Mach number of 0.77. The flutter dynamic pressure for the numerical model is 150.8 psf. The 
uncorrected IS AC generated models indicate flutter occurs between 156 and 163 psf. The flutter frequency of 
the BACT wind tunnel model is approximately 4 hertz. At the same operating condition the flutter 
frequency of the numerical model is 4.16 hertz. The uncorrected IS AC generated models indicate the flutter 
frequency to be approximately 4.22 hertz. 

In terms of the flutter dynamic pressure and frequency at the Mach number for which aerodynamic data 
were available, the numerical model of the BACT system developed herein (and the IS AC-based models) 
gives excellent results. 

Transfer Function Comparisons 

One of the most important measures of model fidelity for control system design is the frequency 
response. In order to effectively design a control system to stabilize a flutter mode the design model must 
accurately characterize the dynamic behavior of the aeroelastic system over a fairly wide range of dynamic 
pressures from stability to neutral stability to instability. 

Figures 5 and 6 show representative frequency responses for the numerical model and the actual wind 
tunnel model. The operating conditions corresponds to a subsonic Mach number and dynamic pressure of 
125 psf, well below flutter. In figure 5 the output is trailing edge accelerometer and the input is trailing 
edge control deflection. In figure 6 the output is trailing edge accelerometer and the input is upper spoiler 
control deflection. The frequency response of an uncorrected IS AC-based model is also presented in figure 5 
for comparison. Note that ISAC cannot model spoilers and so no ISAC comparison can be made in 
figure 6. The experimental transfer function for the upper spoiler was obtained by operating the spoiler 
about a bias. That is, the spoiler was given a static deflection (typically between five and ten degrees) and 
oscillated symmetrically about the bias. The frequency responses of the other accelerometers are comparable 
to those shown here. 

There is excellent correlation between the experimentally obtained frequency responses and those of the 
numerical model. The model clearly captures the key aspects of the dynamic response of the BACT system 
at the subcritical dynamic pressure of 125 psf. There are, however, some discrepancies in the frequency of 
the magnitude peak near 3.5 hertz and in the phase characteristics. However, these differences are typical of 
aeroelastic systems and can be attributed to modeling errors and/or artifacts of using FFT's to generate 
frequency responses from experimental data. 
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Frequency (Hz) 

Figure 5 - Frequency response of trailing edge inboard accelerometer due to trailing edge control: 
q = 125 psf. (solid - experiment, dash - model, dot - IS AC). 




Frequency (Hz) 

Figure 6 - Frequency response of trailing edge inboard accelerometer due to upper spoiler control: 
q = 125 psf (solid - experiment, dash - model). 
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Open-Loop RMS Accelerations 

It is important for the numerical model to accurately characterize the response of the system to 
disturbances because disturbance response determines the control activity required to achieve the desired level 
of closed-loop performance. The disturbance source for the BACT system is wind tunnel turbulence. One 
measure of the degree to which the numerical model characterizes the effects of turbulence is rms 
acceleration. Table 8 presents rms trailing edge inboard acceleration at two dynamic pressures. The data are 
based on normalized dynamic pressure expressed as a fraction of the flutter dynamic pressure, qj. 

Normalization is needed because of the differences in the flutter dynamic pressure from the experiment and 
the numerical model. The reference speed used to scale the turbulence model is 400 fps and is consistent 
with the speed of the flow in the wind tunnel. 


Table 8 - Comparison of RMS Trailing edge Acceleration 


^norm 

RMS Trailing Edge Acceleration (g) 


(psf) 

Experiment 

Model 

% Error 

0.75 q j 

0.0207 

0.0188 

-9.2 

0.90 q f 

0.0340 

0.0350 

2.9 


Note that there is a Mach number mismatch between the experimental data and the model-based data 
since the aerodynamic parameter values in the numerical model are based on data collected at Mach 0.77 and 
the experimental data was obtained for Mach numbers of 0.63 and 0.71 . The good agreement in the response 
level implies that the numerical model can be used to effectively assess rms response. 

Comments 

Based on the accuracy of the flutter properties, the subcritical frequency responses, and rms disturbance 
response, it is reasonable to expect the model to characterize the dynamic response of the BACT system 
over the anticipated range of operating conditions and is appropriate for control system design. However, the 
discrepancies identified above should be considered during control system design and dynamic analysis. For 
example, uncertainty models used in robust multivariable control system design could characterize 
uncertainty in the pitch effectiveness of the trailing edge control surface, the differences in magnitude and 
phase of the system frequency responses, and the differences in flutter dynamic pressure and frequency. 

A MATLAB/SIMULINK* implementation of the wind tunnel model, actuator models, turbulence model, 
and digital controller effects has been developed for the purpose of evaluating and analyzing the dynamic 
behavior of the BACT system. ^ 14 ^ It has been used by several researchers to aid in the design and analysis of 
flutter suppression controllers. A variety of classical, Hoo, p-synthesis, neural network, and adaptive 
controllers have been designed using the numerical model and have been successfully tested in the Langley 
Transonic Dynamics TunnelJ 15 ’ 16,17 ^ 

Concluding Remarks 

The dynamic model of the Benchmark Active Control Technology (BACT) wind tunnel model 
presented herein has many advantages over a purely numerically derived model. It is analytical and 


* MATLAB and SIMULINK are registered trademarks of The MathWorks, Inc. 
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parametric in nature and therefore lends itself to sensitivity and uncertainty analysis. Because the 
aerodynamic effects are represented in derivative form, experimental data can readily be substituted for 
analysis-based data. A major advantage of the modeling approach is that it allows experimental stability and 
control derivative data to be used to model spoiler aerodynamics. The modular form of the model also 
allows various components of the model to be modified or replaced. This is very useful in cases where 
actuator models and turbulence models are modified or updated. 

The BACT model and test data are also being developed as a case study package for educational use. The 
relatively simple structure of the BACT system coupled with the availability of extensive and detailed 
experimental data make the BACT an excellent candidate for additional study of dynamics and control of 
aeroelastic systems. 
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Appendix 


Definition of Symbols for Matrix Form of BACT Equations of Motion 


Vector of generalized coordinates - 



Vector of control surface inputs - 



Vector of disturbance inputs - 



(Al) 


(A2) 


(A3) 


Structural, aerodynamic, and inertial coupling mass matrices - 


M r 


m s he 
s h8 h 


M a S 


-Ci 


cC 


M n 


- f (x)C La 

ci{x)C Ma 




m 

s h6 


(A4) 


Structural and aerodynamic damping matrices - 


2 ^h w h 

0 ■ 

D n = 


0 

2£ e a) e 

a 

cC Ma c((x)C Ma + (c Ma + C Mq ) 


(A5) 


Structural and aerodynamic stiffness matrices - 


K 


S 


s 

K h 

0 ■ 

K n m 

0 

- c ^i 


0 

K °. 


0 



Generalized static aerodynamic force matrices - 

-C, 


Qo - 


c C\t, 


Qt 



(A6) 


(A7) 
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Control effectiveness matrices - 


B o 


-C u -Cf, 1 

&TE 6 US 

E 

\-C L . -C L 1 

&TE 6 US 

m — 

S h6jE Sh &US 

cCm 6je c ^ M dus J 

D \ 

cC M k cC Mi 
°TE °US 


S 66 T e S 8&US 


(A8) 


Numerical Form of Matrices in BACT Equations of Motion 

Structural, aerodynamic, and inertial coupling mass matrices - 


'6.0843 

0.0142' 

M n a 

' 1 .5532 

-0.2712' 

Mg » 

'6.0843' 

0.0142 

2.8017 

a 

-1.3253 

0.2319 

8 

0.0142 


Structural and aerodynamic damping matrices - 


'0.3579 

0.0009' 

D a - 

-4.5840 

1 .0739 ' 

0.0008 

0.1843 

u 

1.490 

-1.7877 


(A9) 


(A10) 


Structural and aerodynamic stiffness matrices - 


'2686 

0 ' 

K n ■ 

0 

-4.5840' 

0 

3000 

a 

0 

1.4900 


Generalized static aerodynamic force matrices - 


Qo e 


O' 

0 


Qt 


-4.584 

1.987 


(All) 


(A 12) 


Control effectiveness matrices - 


B o 


-0.6300 -0.2200' 
-0.0246 0.0573 


B\ 


0 0 
0 0 


-0.0029 -0.00039 
-0.0016 -0.00010 


(A13) 


State Space Matrices for Actuator Models 

The state space form for the actuator models is 


*a ~ A a x a B a & c 


6 

d 

6 


— C a x a + D a S c 


(A14) 


Note that the output of the model consists the acceleration, rate, and position of the control surface. All 
three are used as input to the BACT equations of motion. 
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The state space matrices for the trailing edge control surface actuator are 



-1.8514e + 02 

-2.7324e + 04' 

D _ 

a TE 

'2.7871e + 04' 

A - 

a TE 

' 

0 

0 


-1 ,8514e + 02 

-2.7324e + 04' 


'2.787 le + 04‘ 

c 

TE 

1 

0 

® a TE 

0 


0 

1 


0 


The state space matrices for the upper spoiler actuator are 


-2 ,7880<? + 02 -2 ,6896e + 04' 

ft = 

'3.1 199e + 04' 

« 0 

a us 

0 



'-2.7880e + 02 

-2.6896e + 04' 


'3.1 199e + 04' 

r 

C aus 

1 

0 

D °us 

0 


0 

1 


0 


State Space Matrices for Turbulence Model 

The state space form for the turbulence models is 

x g = \ x g + B g^g 
w 1 

w = Cg x g + D g r) g 


(A 15) 


(A16) 


(A 17) 


Note that the output of the model consists the acceleration and velocity of the downwash in the freestream 
flow. Both are needed as input to the BACT equations of motion. 

The state space matrices for the wind tunnel turbulence model are 


-2.6322e + 02 -1.7321e + 04' 

B* = 

r 

1 0 

g 

0 


(A18) 


-6.1372e + 03 

-7.3026e + 05' 

Do = 

'4.2160e + 0T 

4.2160e + 01 

4.9601e + 03 

g 

0 


The reference speed value for this model is 400 fps because it is close to the experimental test conditions 
and correlates best with experimental data. 
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